7 Be breakup on heavy and light targets 
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In this paper we present all-order quantum mechanical calculations of 7 Be breakup on heavy 
( 208 Pb) and light ( 12 C) targets. We examine the issues concerning the extraction of the astro- 
physical S'-factor 534(0) from the breakup data. We discuss the interplay between Coulomb and 
nuclear breakup, and the importance of higher-order couplings on the cross section. We show that 
nuclear and Coulomb contributions are not separable using the standard angular selection criterion 
as nuclear breakup remains large for small scattering angles, even for the heavy target. However, by 
selecting an upper limit on the relative energy between the final fragments, the contribution from 
the nuclear breakup can be significantly reduced such that Coulomb breakup is the main reaction 
mechanism. We show that the extraction of the asymptotic normalization coefficient may require 
more careful consideration of the nuclear interior than previously used. 

PACS numbers: 25.70.De, 24.10.Eq, 25.60.Gc, 27.20.+n 
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At present, the capture rate 3 He(a,7) 7 Be is now more 
uncertain than 7 Be(p,7) 8 B, both belonging to the pp 
chain and with connections to the solar neutrinos 
Also, the 3 He(a,7) 7 Be reaction is the only important pro- 
duction channel for 7 Li in big-bang nucleosynthesis 0- 
Reference [l| recommends a value of 5*34(0) = 0.53 ±0.05 
keV b. At the low energies of astrophysical relevance, 
capture rates are exceptionally hard to measure directly 
and thus rely on extrapolations. Recent theoretical stud- 
ies claim that the extrapolation of 534(0) can be signifi- 
cantly more uncertain [3|. 

Consequently, measurements with alternative methods 
have been considered. There are two main methods in 
use: i) The Coulomb dissociation method proposes the 
measurement of radiative capture rates from breakup 
data on heavy targets p|, ii) The asymptotic normal- 
ization coefficient (ANC) method proposes the measure- 
ment of the capture rate from transfer reactions [5j or 
breakup reactions Both have been applied success- 
fully to other reactions [|| Q • 

The uncertainties on the direct capture measurements 
for the 3 He(a,7) 7 Be reaction have motivated two exper- 
iments: one at the NSCL, measures the 7 Be breakup 
on 208 Pb at 100 MeV/nucleon, and the other, at the 
Cyclotron Lab at Texas A & M, will use a 7 Be 25 
MeV/nucleon beam on a 12 C target. It is therefore timely 
to consider the reaction mechanisms involved in extract- 
ing 5*34(0) from these reactions. 

The cross sections for measuring the breakup of A — > 
c + x via the Coulomb field of a heavy target T are much 
larger than the low energy direct capture cross sections. 
From the Coulomb dissociation data involving low rela- 
tive breakup energies between the c + x fragments, one 
can extract the inverse reaction c + x — > A at astro- 
physical energies 0, Q . There are three main breakup 
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mechanisms which can complicate this relation: i) nu- 
clear breakup, present whenever the projectile gets close 
to the target, ii) E2 transitions, may be negligible in ra- 
diative capture, but are typically significant in Coulomb 
dissociation 9], iii) final state interactions (continuum- 
continuum (CC) couplings) that distort the final energy 
spectrum of the emitted fragments [T^j ■ All of these en- 
tangle the information on the capture reaction. In this 
paper we examine these various issues for the particular 
case of the 7 Be breakup, making use of the continuum 
discretized coupled channels (CDCC) method, reviewed 
in Ref. Within CDCC, nuclear and Coulomb are 

consistently included and the contribution from CC cou- 
plings can be explicitly explored. 

There have been several applications of the Coulomb 
dissociation method to other astrophysically important 
breakup reactions. The nuclear and Coulomb interplay 
has been studied for the 8 B — > 7 Be+p breakup reac- 
tion 0, 0, ITil ITU . Here, Coulomb dissociation dom- 
inates for small scattering angles, which can be related 
semi-classically to large impact parameters outside the 
range of the nuclear force. Studies of 7 Li and 6 Li res- 
onant breakup have suggested that for these nuclei, at 
high energies, the Coulomb and nuclear breakup can be 
well separated into angular regions, and the nuclear ef- 
fects are small at higher energies and forward angles . 
Conversely, earlier studies of 7 Li breakup have shown 
that experimental data differ significantly from Coulomb- 
only calculations at forward angles which indicates that 
the nuclear forces are important even at forward angles 
[T7L flgl TTg| . In addition, coupled channels calculations 
of sequential breakup via the 3 + resonance in 6 Li have 
been performed, and again nuclear effects were found to 
be important even at forward angles 20] . 

From these studies of similar reactions it is not possi- 
ble to derive the importance of nuclear effects for the 7 Be 
breakup reactions of interest here. We therefore perform 
fully-quantum mechanical calculations, in the CDCC 
framework using the coupled-channels code FRESCO [2lJ , 
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of 7 Be breakup on 208 Pb and 12 C targets at energies of 
100 and 25 MeV/nucleon respectively. 

The breakup of 7 Be into a+ 3 He by the interaction 
with a target consists of a three-body (two-body pro- 
jectile + target) problem. In the CDCC method, the 
breakup of the 7 Be projectile is treated as an excitation 
into the a+ 3 He continuum, discretized into J\f bins. For 
7 Be, with two bound states below the a+ 3 He breakup 
threshold, this requires a J\f + 2 coupled channels prob- 
lem. The physical inputs are then: an a+ 3 He potential 
which binds the 7 Be, and for each cluster, an interaction 
with the target, which can be complex to account for loss 
of flux to unaccounted channels. 

The a+ 3 He potential used here is that of Buck et al. 
[22|. This potential consists of a central and spin or- 
bit terms of Gaussian form. The parameters were fixed 
using the binding energies of the bound states and the 
positions of the resonances for 7 Be and 7 Li, the 7 Li 
charge radius, quadrupole and octupole moments, and 
B(E2:3/2" -> 1/2"). We increased slightly the spin- 
orbit potential depth to fit exactly the 7 Be binding en- 
ergies for the 2p 3 / 2 ground state (S a =1.587 MeV) and 
2pi/ 2 first excited state (5 Q =1.158 MeV), since these fix 
the correct asymptotics of the bound state wavefuntions. 
The bound states have a node because the lower orbitals 
are excluded due to Pauli blocking. This potential pro- 
duced the resonances in the /-waves at approximately 
the correct energies. In Ref. [2^, the potential depth 
was adjusted to fit the s-wave phase shifts, which also 
gave small d-wave phase shifts as required. Therefore, 
we use a parity dependent potential, with a depth of 
V dd = 83.77 MeV for the p- and /-waves, and Kven = 
66.10 MeV for the s- and d- waves. The spin orbit depth 
is Vis = 3.8 MeV, and the range for all three potentials 
is i? dd = -Rcvcn = Ris = 2.52 fm. The Coulomb radius 
is i? c =3.095 fm. 

The calculated B(E2:3/2~ -> 1/2 - ) for this interac- 
tion is 19.0 e 2 fm 4 , consistent with Ref. j2^|. The Buck 
model produces a rms charge radius for 7 Be of 2.62 fm, 
in reasonable agreement with the experimentally deter- 
mined charge radius of 2.52 ± 0.03 fm j24[. 

Optical potentials are important ingredients in our cal- 
culations. The cluster-target potentials are fixed by elas- 
tic scattering data at the energy carried by that cluster 
in the projectile. We require 3 He+ 208 Pb at 300 MeV 
and a+ 208 Pb at 400 MeV. The nearest available energy 
for elastic data was at 217 MeV for 3 He+ 208 Pb [H] and 
340 MeV for a+ 208 Pb 26j. We require 3 He+ 12 C at 75 
MeV and a+ 12 C at 100 MeV. We use the potential for 
3 He+ 12 C at 72 MeV from Ref. [27| but neglect the spin- 
orbit term, and the potential for a+ 12 C at 104 MeV from 
Ref. [2^| . Different sets of potential parameters for the 
fits to the elastic data are reported in the literature. The 
sensitivity of the breakup cross section to the various po- 
tentials was found to not significantly effect the breakup 
cross section (less than 2% at 0° for both targets). 

In the present calculation, the 3 He+a continuum in- 
cludes partial waves for £ < 3 and each j 7r discretized 
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FIG. 1: J-distribution of the cross section for the Pb tar- 
get. The total breakup cross section is shown by the solid line 
while the broken lines give the nuclear (dotted) and Coulomb 
(dashed) contributions. The top scale relates the angular mo- 
mentum to the impact parameter via the semi-classical rela- 
tion J = Kb, where K — 15 fm -1 . The insert is the same 
plot on a linear scale. The arrows are discussed in the text. 
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FIG. 2: J-distribution of the cross section for the C target. 
The lines have the same meaning as Fig. Q and the impact 
parameter scale on the top axis uses K = 4.8 fm -1 . 



into 10 bins up to 20 MeV, equally spaced in momen- 
tum. Due to the sharp resonance, the fa 12 partial wave 
has a resonance bin with the width of the resonance, and 
two bins below the resonance bin. Above the resonance, 
we include 8 bins up to 20 MeV for the lead target and 9 
bins up to 30 MeV for the carbon target. Convergence of 
this discretization of the continuum was checked by dou- 
bling the number of bins, which only affected the total 
breakup cross section by 2%. The continuum bins are 
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integrated out to a maximum radius rbi n = 50 fm. Po- 
tential multipolcs arc included for q < 2 as octupole tran- 
sitions are negligible. Using interpolation, partial waves 
up to L — 14000?i are included for the lead target and 
up to L — 2000?i for the carbon target. The coupled 
equations are solved using i? max = 1000 fm. We use non- 
relativistic kinematics which, at 100 MeV/nucleon on the 
lead target, will introduce a 3% error in the momentum. 

In Fig. H (Fig- 0! we show the J-distribution of the 
breakup cross section for the 208 Pb ( 12 C) target. The im- 
pact parameters that correspond to each partial wave, us- 
ing the semi-classical relation J = Kb, are shown across 
the top scale. The total breakup cross section (solid line) 
is shown along with the nuclear breakup (dotted) and 
Coulomb breakup (dashed). The sum of the radii for 
the projectile and target is marked on the figures by the 
down-pointing arrow. 

Expectedly, we see that for both light and heavy tar- 
gets, the nuclear breakup is the dominant process for the 
lower partial waves and Coulomb breakup dominates the 
higher partial waves. The Coulomb breakup, suppressed 
at low partial waves due to the nuclear absorption, peaks 
around the sum of the projectile and target radii, then 
falls off slowly with increasing J. The nuclear breakup 
is small for the low partial waves since the imaginary 
part of the nuclear potential is removing flux from elas- 
tic breakup to other channels. It peaks sharply around 
impact parameters corresponding to surface collisions of 
the two nuclei, then falls off rapidly. 

To extract the ANC from which 534(0) can be deter- 
mined, breakup data from a range of targets can be used. 
The fundamental requirement is peripherality. A simple 
sum of radii would imply that there should be no con- 
tribution to the breakup from impact parameters below 
9.8 (5.3) fm for the lead (carbon) target. The plots of 
Figs. Hand |21 suggest that this may not be the case. The 
breakup cross section for impact parameters below the 
sum of radii is 28% (16%) of the total breakup cross sec- 
tion for the lead (carbon) target. 

The slower fall off for Coulomb breakup means that 
the Coulomb dominates beyond a certain partial wave, 
where the nuclear breakup is negligible and interference 
effects are small. The impact parameter, beyond which 
nuclear effects can be considered negligible (up-pointing 
arrow on figures), is around 16 (13) fm the lead (carbon) 
target, values much larger than the sum of the projectile 
and target radii. In the Coulomb dissociation method, 
data are typically taken at forward angles since then 
it is assumed to be nuclear free. For pure Rutherford 
trajectories, there is a relationship between the impact 
parameter and the scattering angle: for 7 Be+ 208 Pb at 
100 MeV/nucleon, an impact parameter of 16 fm corre- 
sponds to a center-of-mass scattering angle of 2.5°; for 
7 Be+ 12 C at 25 MeV/nucleon, 13 fm corresponds to an 
angle of 1.4°. However, the determination of these cut- 
off angles is rather simplistic. One should note that the 
Coulomb breakup cross section here differs significantly 
from the semi-classical Coulomb dissociation cross sec- 
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FIG. 3: Angular distribution of cross sections for 7 Be elastic 
breakup on 208 Pb at 100 MeV/nucleon. The CDCC (solid) 
calculation is broken down into nuclear (dotted) and Coulomb 
(dashed) contributions. The sum (dot-dashed) represents the 
incoherent sum of the nuclear and Coulomb contribtutions, 
while the solid line is the coherent sum. 
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FIG. 4: Angular distribution of cross sections for 7 Be elastic 
breakup on 12 C at 25 MeV/nucleon. The lines have the same 
meaning as Fig. [3] 



tion [2^|. This is mainly due to the finite size of the 
projectile which has been previously pointed out in the 
CDCC calculations of Ref. [l^. 

We can directly examine the angular distribution of the 
breakup cross section on the lead (carbon) target, shown 
in Fig. (Fig. 0J. The nuclear (dotted) and Coulomb 
(dashed) contributions to breakup are plotted along with 
the coherent sum (solid). The interference between nu- 
clear and Coulomb breakup is shown by the difference 
between the incoherent sum (dot-dashed) and the coher- 
ent sum (solid). We see that for both targets the nuclear 
breakup is the dominant process, having a diffractive na- 
ture which peaks at zero degrees. Contrary to expec- 
tations, we see that the nuclear and Coulomb breakup 
cannot be separated into angular regions and both have 
to be considered. 
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FIG. 5: Angular distribution of the breakup cross section in 
Fig. [3] which includes only the lowest two energy bins from 
each f set. This gives the relative energy between the a and 
3 He fragments an upper limit of approximately 1 MeV. 
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FIG. 7: Angular distribution of the cross section for 7 Be 
breakup on 208 Pb at 100 MeV/nucleon. The different calcula- 
tions of the cross section are: CDCC (solid), DWBA (dashed), 
and a CDCC calculation without continuum-continuum (CC) 
couplings. 
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FIG. 6: Angular distribution of the breakup cross section 
showing the dipole and quadrupole contributions. Both nu- 
clear and Coulomb are included. As with Fig. this includes 
only the lowest two energy bins from each set. This gives 
the relative energy between the a and 3 He fragments an upper 
limit of approximately 1 MeV. 



However for the heavy target, the nuclear breakup con- 
tribution can be reduced by imposing an upper cut on 
the relative energy between the a and 3 He fragments. To 
show the effect of this energy cut on the angular distribu- 
tion we sum up the angular cross section from the lowest 
two energy bins from each j n set. This restricts the max- 
imum final state relative energy to approximately 1 MeV 
(0.846 MeV for all j* except the /V/2 which has a maxi- 
mum energy of 1.268 MeV). The angular distribution of 
the breakup cross section with this energy cut is shown 
in Fig. [3] We now see that Coulomb breakup is dom- 
inant except for the extreme forward angles (less than 
0.5°). To give an estimate of the E2 contribution to the 
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FIG. 8: Angular distribution of the cross section for 7 Be 
breakup on 12 C at 25 MeV/nucleon. The lines have the same 
meaning as Fig. |7| 



breakup cross section, we plot the dipole and quadrupole 
contributions to the breakup (solid line in Fig. which 
includes nuclear and Coulomb breakup) . Once again the 
maximum relative energy has been restricted to approx- 
imately 1 MeV so there is a Coulomb dominated region 
up to # c .m. = 3°, as shown in Fig.^l We see that in this 
region, dipole transitions dominate. 

For lighter targets, the nuclear contribution to breakup 
is usually large and Coulomb breakup is considered small. 
We see from Fig. |J| that this is the case above 6°. Below 
6° the coherent sum of the nuclear+Coulomb breakup 
(solid) is very similar to that of nuclear (dotted) breakup 
alone, although this does not mean that Coulomb ef- 
fects are negligible. Below 6°, the nuclear and Coulomb 
breakup contributions interfere, through strong couplings 
in the continuum, to give a coherent sum which is similar 
to that of nuclear alone. 
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The effect of couplings on the angular distribution is 
shown in Fig. [7| (Fig. HJ for the lead (carbon) target. 
The solid line is the full CDCC calculation which in- 
cludes all couplings to all orders. The DWBA calcu- 
lation (dashed) only includes first order couplings be- 
tween ground and excited states. The dotted line shows 
a subset of the CDCC calculation which neglects the 
continuum-continuum (CC) couplings. We see that for 
both targets the couplings reduce the cross section sig- 
nificantly, less so for the heavier target, where the exper- 
iment was performed at a higher beam energy. When ex- 
tracting astrophysical quantities from this breakup data, 
it is important to go beyond first-order reaction theories. 

In conclusion, wc have shown that for these breakup 
measurements, contributions from the interior are sig- 
nificant. Selected cuts on the data need to be consid- 
ered for the extraction of the ANC, which requires a 
peripheral collision. In addition, the relationship be- 
tween the ANC and the astrophysical 5-factor relies on 
a first order DWBA approach, but we have shown that 



the DWBA cross section dramatically overestimates the 
breakup cross section. It is not clear whether this method 
will be helpful in improving 5*34(0). We have shown that 
by performing 7 Be breakup experiments on heavy tar- 
gets, it is not possible to completely eliminate nuclear ef- 
fects through selecting 6* c . m . smaller than a critical value. 
However, by selecting an upper limit on the relative en- 
ergy between the final fragments, the contribution from 
the nuclear breakup can be significantly reduced such 
that El Coulomb breakup is the main reaction mecha- 
nism. The extraction of 5*34(0) using the Coulomb disso- 
ciation method may be promising, but as yet, the uncer- 
tainties due to the nuclear and E2 contributions, along 
with the experimental uncertainties in the measurement, 
makes it unclear whether this method will be able to im- 
prove on the radiative capture measurements. 
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